A Numerical Analysis to the n and K 
Coupled— Channel Scalar Form-factor 

Wei Liu, Hanqing Zheng and Xiao-Lin Chen 

1) Department of Physics, Peking University, Beijing 100871, P. R. China 

Abstract 

A numerical analysis to the scalar form-factor in the tttt and KK 
coupled-channel system is made by solving the coupled-channel dis- 
persive integral equations, using the iteration method. The solutions 
are found not unique. Physical application to the tttt central produc- 
tion in the pp — > ppirir process is discussed based upon the numerical 
solutions we found. 
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The I=J=0 channel tttt interactions are of great physical interests. Be- 
cause the interaction of the I=J=0 channel is very strong, the input bare 
singularities for a given model can be severely renormalized and distorted 
by the strong attractive force and extra dynamical singularities may be gen- 
erated 0. Therefore the I=J=0 channel affords an ideal test ground for 
models of strong interactions. Also, the lightest glueball is expected to lie in 
the I=J=0 channel, as predicted by the lattice QCD calculations. Therefore 
a detailed study to the dynamics in this channel becomes especially inter- 
esting. The low energy I=J=0 tttt system manifests itself in various produc- 
tion processes intensively measured by experiments, i.e., from ttN — > ttttN, 
77 — ► tttt to J/^ — ► </>7T7r, etc., with the center of mass energy, \/s, ranging 
from the tttt threshold to a few GeV. When y/s exceeds the KK threshold, a 
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single-channel analysis to the tttt system becomes inadequate, and instead, a 
coupled-channel analysis of tttt and KK system has to be made. Among var- 
ious tttt production processes particularly interesting cases are those I=J=0 
final states which are generated weakly. In such circumstances, the final 
state particles will not scatter back to the initial states, and therefore the 
complicated dynamics involved is considerably simplified without loss of in- 
formation on the I=J=0 final state interactions. For example, in such a case 
the Watson-Migdal's theorem on final state interactions applies. The tttt fi- 
nal states "weakly" produced can be again categorized into two classes: One 
is that the tttt production vertex contains left-hand singularities like in the 
case 77 — > tttt [] and the another is not, like in the case K — > tttt and in the tttt 
central production process pp — > ppmr. Physical situation in the absence of 
left hand singularities is further simplified, since the dynamical complexity 
from the production vertex is removed which would otherwise disturb our 
analysis on the tttt (final state) interaction itself. In such a simplified situa- 
tion, it then becomes reasonable to assume that the tttt production amplitude 
is factorized as a product of a form-factor-like quantity, which we denote as 
A, containing all the dynamical singularities from the right-hand cut from 
tttt final state interactions, and the production vertex which is a smooth an- 
alytic function of s on the complex s-plane except possibly at infinity since 
it contains neither the left-hand nor the right-hand singularities. 

For a coupled-channel system of tttt and KK, the spectral representation 
of the form-factor, A = (Ai,A 2 ), satisfies the following relation, 

ImAi = A 1 p 1 T{ r 1 + A 2 p 2 T 2 \ , m 
lmA 2 = A 2 p 2 T+ + A lPl T+ , 1 > 

whereas the unitarity relation of the scattering matrix T reads, 

ImT = TpT+ , (2) 

where p = diag(p\, p 2 ) is the matrix of the kinetic phase-space factor. The T 
matrix may contain left-hand cut but A does not. Especially, A is analytic 
on the entire physical sheet of the complex s plane except on the cut along 
the real positive axis starting from 2tt threshold. The form-factor A has the 
same analytic structure as the scalar form-factor Q and the two are different 

A recent publication on related subject can be found in Ref. @. 
2 For the definition of the scalar form- factor, see for example Ref. M. 
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only up to a polynomial. Eqs. ([[]) and (0) are assumed to be correct down to 
the lowest threshold, in the absence of anomalous thresholds. That means 
when 4m^ < s < 4m\ Eq. (JTJ) takes the form, 

\mA l = A 1 p 1 T^ 1 , 

lmA 2 = A lPl T+. [6> 

From previous discussions, it is realized that the physical problem of 
studying final state interactions in the production process without left-hand 
singularities is reduced to the mathematical problem of solving Eq. (|I|) (and 
Eq. (|3D), provided that the T matrix is known. In the single channel case 
the analytic solution of the form factor can be obtained. The spectral rep- 
resentation of the form-factor is given by the first equation in Eq. (|3|) from 
which the classical Omnes solution can be established, 

A(s) = P(s) exp (£ r ^ S>) , da) , (4) 

where P(s) is a polynomial and 8^ is the tttc scattering phase. The Omnes 
solution is remarkable in relating the form-factor to 5 n . The Eq. (^) has 
been used in Ref. || to study the final state interactions in K — > 2tt and 
pp — > ppiTTc systems. In the coupled-channel case, however, no analytic 
solutions of Eq. (|I]) can be obtained in general. In the following we will 
instead study the coupled-channel system by numerical method ||. 

Since the coupled-channel form-factor A is analytic on the physical sheet 
of the complex s plane except for the cut along the real axis, assuming the 
T matrix is known we can search for solutions of the amplitude A in Eq. ([I]) 
by solving the following dispersion relation: 

A= i , A^ry) 

7r jr s' — s — ie 

where the integration is performed on the unitarity cut R, starting from 4m^ 
to oo. The Eq. (|5|), according to Muskhelishvili ||, contains two fundamental 
solutions, 0i and 02- Assuming that n behaves as n — > s~ Xn as s — > oo 
one has 

y^Xn = —[argdetS] R = -(<J ff (oo) + S K (oo)) , (6) 
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where S is the coupled-channel S matrix. According to ||, any solution of 
the integral equation (|5]) can be written as a linear composition of the two 
fundamental solutions, 

= ^P„( S )0 n , n=l,2 (7) 

n 

where P n (s) are polynomials of s. The polynomials are not determined from 
analyticity alone. Other physical input has to be implemented to fix their 
coefficients. For example, chiral perturbation theory can afford an expansion 
of the scalar form-factor in powers of s when s is small 0. Since there 
can be many solutions of Eq. (^j) and we notice that, since the numerical 
integration in Eq. (j5|) has to be truncated somewhere (denoted as A below), 
all the information on the asymptotic behaviour are lost. As a consequence, 
it is difficult to distinguish the so-called fundamental solutions from others, 
since the difference between the two essentially comes from their asymptotic 
behaviours. Therefore in the present numerical scheme instead of searching 
for the fundamental solutions, we follow the recipe of Ref. [Q, that is to search 
for two linearly independent solutions, A 1 and A 2 which are normalized at 
s = as, 

A{(0) = 1, 4(0) = 0, (8) 

and 

4(0) =0, 4(0) = 1. (9) 

In the following we use the im and KK coupled-channel fit of the T 
matrix from Au, Morgen and Penington [[/j] as an educative example to solve 
the coupled-channel dispersive integral equation (|5]). The integral in Eq. (|J) 
is truncated at A ~ l.bGeV . We find that the influence of the cutoff is rather 
local, i.e., it has little effects on the behaviour of the solution in a large region 
of s, even close to A. One can also use a mild regulation function instead of 
the truncation at A and the result is essentially the same except at s ~ A. 
We use the iteration method to solve Eq. (|J). That is, 

, s M Real fA( n )(s')pT+(s / ) 
A^\s) = A(0) + £ / 1 } ± M , (10) 

7T J4m| S'{S' — S — It) 

where 

A(0) = (1,0) or (0,1). (11) 
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A once-subtracted form of the dispersive integral in the above Eq. ( |T0D is 
helpful in incorporating the boundary conditions, Eqs. @ and @. The 
routine converges rather rapidly (after about 20 steps), which confirms the 
claim in Ref. H and the solutions are not unique due to the reason already 
mentioned above. This means that the iteration may converge to different 
solutions depending on different initial values for the iteration. In Fig. p] and 
Fig. |^ a few examples generated from our numerical recipe are shown. We see 
from Fig. [I] that all the phases are identical to the phase of T\\ below the KK 
threshold as required by the final state theorem. Above the KK threshold, 
however, the phases from different solutions of A\ can be very different and 
also deviate from the phase of T\\. In Fig. || the corresponding magnitude 
of the different solutions is also shown. We see that around lGeV region 
there can be zeros or dips to compensate the peak generated by / (980). The 
physical discussion on the necessity to introduce the protective zero appearing 
in the nn production processes can be found in Ref. H [|. In the present 
numerical approach there is generally no difficulty to pick suitable solutions, 
from all, to fit experimental data. An example is shown in Fig. |3] where the 
7T7T production cross-section in the pp — > pprnr process is reproduced. The 



cross-section can be expressed by [1C 



da (.-4n£)W 



where F(s) is a linear combination of solutions A\ and A\, 

F(s) = a 1 A\ + a 2 A\ . (13) 

It is assumed here that ttti produced in the pp — > ppim process are generated 
from the fusion of Pomeron pairs and no left-hand singularity is involved in 
the production vertex. Therefore a% and ot 2 (which contain information about 
the 7T7T production vertex) appearing in the above equation are expected to 
be polynomials with a weak dependence on s. 

In fact, the T matrix given in Ref. J7[ is obtained by a multi-pole K 
matrix fit in which the background effects including the left-hand cut one 
are simulated by a polynomial. Therefore, the T matrix does not contain any 
left-hand singularity. Under such an approximation it is not really necessary 
to search for the numerical solutions of A\ in order to make use of Eq. (|13|) . 
Since A\ has to be a linear combination of Tu and T 2 \ when neglecting the 
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left-hand singularities in T, one can use instead of Eq. ([13]) the following 
form, 

F(s) = apu + a' 2 T 21 , (14) 

to fit the experimental data, as originally done in Ref. [[F]. The method pre- 
sented in the current note applies to more general cases when the fine details 
of the left-hand singularities of the scattering T matrix are carefully taken 
into account. Also our program can be easily extended to the situation when 
the left-hand cut of the production vertex is included. We will investigate 
these more realistic cases in future. 
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Figure 1: The solid, dashed and dot-dashed lines represent the phases of 
different solutions of the A\ amplitudes, the corresponding magnitude of the 
amplitudes are depicted in Fig. |[ The dotted line represents the phase of 
Tu from the K3 fit of Ref. M. 
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Figure 2: The magnitude of the different solutions of the A\ amplitude. 
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Figure 3: The effective mass distribution of pion pairs in pp — ► ppinr vs. 
M = y/s. The data are from [||. 
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